Quasilinear spin voltage profiles in spin thermoelectrics 
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Recent experiments show that spin thermoelectrics is a promising approach to generate spin 
voltages. While spin chemical potentials are often limited to a surface layer of the order of the spin 
diffusion length, we show that thermoelectrically induced spin chemical potentials can extend much 
further in itinerant ferromagnets with paramagnetic impurities. In some cases, conservation laws, 
e.g., for a combination of spin and heat currents, give rise to a linear spin voltage profile. More 
generally, we find quasilinear profiles involving a spin thermoelectric length scale which far exceeds 
the spin diffusion length. 



Introduction. — The field of spintronics 0, Q is fueled 
by the vision of more efficient microelectronic devices 
which exploit the spin degree of freedom a of electrons 
in addition to their charge e. To date, attention focused 
largely on electric phenomena. Recently, much considera- 
tion is also given to spin-dependent thermoelectric effects 
which opened up the new research field of spin caloritron- 
ics 0] . A potential milestone of this emerging field were 
the recent observations of the spin Seebeck effect in ferro- 
magnets , where a spin voltage develops in response 
to an applied temperature gradient. This could be valu- 
able for spintronics applications as the operating princi- 
ple of a spin battery. Indeed thermally driven spin injec- 
tion has recently been realized experimentally [7|. 

Besides such practical implications, the spin Seebeck 
effect also raises theoretical puzzles. Remarkably, the 
spin chemical potential is observed to vary (quasi)linearly 
along the sample jj-Q. This is surprising since unlike 
charge, the electronic spin is generally not conserved due 
to spin-orbit interactions or magnetic impurities. While 
charge conservation implies that the chemical potential 
drops linearly, spin non-conservation would seem to im- 
ply that the spin chemical potential should vary expo- 
nentially on the scale of the spin diffusion length ||. 

Proposed mechanisms for the linear spin-voltage pro- 
file are based on spin pumping between ferromagnet and 
contacts 0- 11 1 , with previous theoretical works |9rll~2| fo- 
cusing on ferromagnetic insulators or magnon drag. In 
this paper, we consider the contribution of the conduc- 
tion electrons to the spin Seebeck effect in a ferromag- 
netic metal. In particular, we discuss the question under 
what circumstances a (quasi)lincar variation of the spin 
chemical potential can be obtained. Although it is well 
known that elastic spin-flip scattering leads to an ex- 
ponential decay of the spin chemical potential near the 
edges of the sample Q , we find that the situation changes 
qualitatively when additional inelastic spin-flip scattering 
is present as, e.g., arising from the exchange interaction 
with magnetic impurities. 

As a model system, we consider inelastic spin-flip scat- 
tering effected by the exchange interaction with partially 
spin-polarized (paramagnetic) impurities. In the absence 
of elastic spin relaxation of the conduction electrons or 



without additional relaxation mechanisms for the impu- 
rity spins, a truly linear dependence of the spin voltage 
can be recovered due to the interplay of a conservation 
law for a combination of spin and heat currents with the 
spin asymmetry inherent to inelastic spin-flip scattering. 
Remarkably, even in the more general case with addi- 
tional direct spin relaxation for both conduction electrons 
and impurities, we find a quasilinear variation of the spin 
voltage whose characteristic scale is much larger than the 
spin diffusion length. 

Another question raised by the experiment is the origin 
of the observed spin Seebeck effect |4|] . In the absence of 
response functions which are off-diagonal in spin, linear 
response relates the spin-resolved charge and heat cur- 
rents, J u and Jq ct , to the gradients in the spin-resolved 
chemical potentials and temperatures, \x a and T a , 
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Here, o~ a denotes the (spin-resolved) conductivity, S a the 
thcrmopower, and n a the heat conductivity. It is natural 
to ask whether the observed spin Seebeck effect indeed 
originates from the spin thcrmopower S- = — Sj,. We 
find that this occurs in special limits but is generically 
not the case. 

Model. — We consider an itinerant ferromagnet where 
the exchange interaction polarizes the electron spins 
through the mean field h and the electron dispersion 
is spin split, e pa = e p — ah. Besides elastic potential 
scattering we assume that the electrons are subject to 
inelastic spin-flip scattering (with rate T S J) on partially 
polarized magnetic impurities. The electron spin can also 
relax directly (I^ ir ), e.g., via spin-orbit interaction. Sim- 
ilarly, the impurities can relax their spin directly (1^ r ), 
e.g., via the combined action of spin-orbit interaction and 
phonons [fij ] in addition to spin-flip scattering with the 
conduction electrons (I^ f ). 

Our discussion of the spin thermoelectric properties of 
the ferromagnet is based on the kinetic equation 
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for the electronic distribution function n prj . Here, v pcr = 
V p e p(T , and F is an externally applied force. Elastic 
scattering and inelastic spin-flip scattering are described 
through the collision integrals 

S P<r = J2 W PP' 6 ^ - £ P'-') (V'-V) 
p'a' 

p' 
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where g a ± = g+(l±a)/2. Here, W aa results from clastic 
impurity scattering and can be spin dependent as, e.g., 
for spin polarized impurities, W" denotes the elastic di- 
rect spin relaxation, and W st denotes the strength of the 
inelastic spin-flip scattering with paramagnetic impuri- 
ties. For definiteness, we consider isotropic scattering 
due to pointlike scatterers. Since paramagnetic impuri- 
ties in a ferromagnet will be at least partially spin polar- 
ized, the impurity -levels split by an energy of e^. We 
treat the impurity spin as classical, i.e., we describe devi- 
ations from the maximal spin orientation with a bosonic 
distribution function g, which is governed by 
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W°± denotes direct relaxation of the impurity spin. 

Transport coefficients. — The transport coefficients a a , 
Scr, and K a in Eq. ([T]) can be determined as usual by eval- 
uating the spin-resolved charge and heat currents from 
the Boltzmann equation ((2|). Multiplying Eq. ((2]) by ve- 
locity v pcr and v pcr (e pcr — //o) (where P o is the equilib- 
rium chemical potential) and summing over momenta, 
we find a a = e 2 N a D a for the conductivity. Here, N a 
is the electronic density of states and D a = v^/aT® 
is the diffusion coefficient in terms of v a = v PF ^ a and 
r° = r£ + r^ f + r^ ir , where r^ r is the transport rate for 
spin-conserving elastic scattering, T^ lr = N s { Po )W aa , 
and F^ f = N s (fiQ + cre.;)(W" 3f /smh -jrpp) the rates for elas- 
tic and inelastic spin-flip scattering, respectively. The 
thermopower and the thermal conductivity obey the 
Mott relation S a = (iT 2 k 2 B T /3e)(d\xio- a { P o) / d^o) and 
the Wiedemann- Franz law K a = (Tr 2 kgT/3e 2 )a a . 

Continuity equations. — Information about the spatial 
dependence of (spin) chemical potential and (spin) tem- 
perature is contained in the continuity equations for 
the spin-resolved charge and heat currents. These can 
be found in the usual manner by a momentum sum 
over the Boltzmann equation (|2|) as well as over Eq. 
d2]) multiplied by (e pa — /i ). We evaluate the result- 
ing expressions by linearizing the collision integrals using 
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Eliminating the impurity distribution function g using 
Eq. ((J]), we find, in the static limit for the charge cur- 



heat current Jq = Jq-j- + Jqj., and the spin heat current 
J<2spin = Jgt _ J Qi' tnc continuity equations 
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Here we defined the combined scattering rates T a = 
T S J + Tf r and f CT = oFf + Tf T as well as the coeffi- 
cients a = rf /(rf + rf ) and a Q = r*/(r? + rf ). 

The coefficient a measures the rate of direct relaxation of 
the impurity spins Tf 11 ' = Wf" — Wf" relative to spin-flip 
scattering with the conduction electrons, 
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i.e., a varies between zero when the impurity spin relaxes 
only via its interaction with the conduction electrons and 
unity when direct relaxation dominates. Furthermore, 
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are the electronic densities and currents. Note that the 
continuity equations contain heat densities or densities of 
even higher order in energy such as P qi due to the energy 
dependence of the spin- flip scattering rates T S J and r^ lr . 

The continuity equations (J5j show that the charge cur- 
rent is always conserved while the heat current is con- 
served only in the absence of inelastic scattering (e; = 0) 
or of direct relaxation of the impurity spins (a = 0). 
Both inelastic spin-flip scattering and direct spin relax- 
ation act as sources or sinks for the spin currents. 

It is instructive to consider first the spin chemical po- 
tential in the absence of thermoelectric effects. Then, 
we can focus on the continuity equations for the charge 
and spin currents, and neglect the energy derivatives T' a 
of the scattering rates. Inserting Eq. ([1]) with S a = 0, 
expanding Pc = N a ([A a — P o), and rewriting the equa- 
tions in terms of the charge and spin chemical potentials 
P± = (V 2 )[0"t - Po) ± (Ml - Mo)] yields the coupled 
diffusion equations 
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Here, wc defined a ± = (a 1 -±a l ) and T%_ = A^r^iTV^Pj*. 
Note that T_ is nonzero only due to the inelastic na- 
ture of spin-flip scattering, as becomes evident by the 
more explicit expression T s _l = \Ni-(uo)Ni(ny + ti) — 
N^N^o-e^iW^/smh-^) Q. For Eq. © one 
finds a linear and an exponential solution. Determining 
the spin voltage profile for a fcrromagnct of length L with 
a voltage Afi applied along the x-direction (correspond- 
ing to the boundary conditions fx+(x = ±L/2) = ±A/i/2 
and fj,-(x = ±L/2) = 0) one finds for a weak ferromag- 
net Ci,h <C ep at small temperatures fc^T <C ep in the 
bulk of the sample: /i_ ~ — ^ + oT^/r+. This implies, 
that the linear spatial variation of the charge chemical 
potential fi + due to charge conservation is necessarily 
accompanied by a linear variation of the spin chemical 
potential when is nonzero. Interestingly, the linear 
contribution to /i_ is entirely determined by the collision 
integral and thus insensitive to the response coefficients 
such as the conductivity. 

Spin thermoelectrics. — This mechanism does not carry 
over directly to spin thermoelectric effects since the latter 
are generically investigated in the absence of charge cur- 
rents. Similarly, an ideal spin battery would also be oper- 
ated without charge currents. Although charge conserva- 
tion becomes inconsequential under these circumstances, 
we find that a linear or quasilinear spatial dependence 
of the spin chemical potential due to inelastic spin-flip 
scattering persists under quite general circumstances. 

For reference, we first review what happens without 
inelastic spin- flip scattering (r sf — 0), as considered in 
Q. Then, both the spin chemical potential and the spin 
temperature T_ = (T^ — TjJ/2 exhibit a purely expo- 
nential dependence, despite the presence of the conser- 
vation law V • Jq = in addition to V • J = 0. The 
reason is that in the equations for the linear modes, 
/j_ and T_ decouple completely from /i+ and T+ = 
(l/2)[(T t - T ) + (Tj. - T )]/2. Indeed, expanding the re- 
maining continuity equations ([5]) around fj-o and Tq yields 
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where 7 = (7r 2 /3)/cgXo, implying that the only sinks and 
sources of the (heat) spin current are [i- and TL. 

As in the absence of thermoelectric effects, spin-flip 
scattering will mix the equations for the charge and spin 
components of chemical potential and temperature. This 
allows for the thermoelectric generation of a linear spin- 
voltage profile in several circumstances. Expanding the 
continuity equations (|5|) around [Lq and To as before, we 
find, in addition to V • J = 0, 

V • J spin = -2(f_/x + + f + 7 f'_T + + 7 f ; + T_) (10) 
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Importantly, jj, + and T + now also act as sinks and sources 
for the currents when F_ is nonzero. Consequently, the 
equations for /i + and T + no longer decouple from those 
for /i_ and T_ and the (quasi)lincar modes involve also 
/j_ and T_. 

A strictly linear spin voltage profile emerges when a 
second current is conserved in addition to the charge cur- 
rent J. If direct spin relaxation is negligible for both 
impurity spins and conduction electrons, the mixing of 
charge and spin components is accompanied by conser- 
vation of the heat and spin currents, V- Jq = as well as 
V • J S pin = leading to a linear spin voltage in the bulk of 
the sample (see discussion below). More surprisingly, a 
linear spin- voltage profile even persists in the presence of 
direct spin relaxation of the impurity spins. At first sight, 
this seems unlikely since both the heat current Jq and 
the spin current J sp i n are no longer conserved. However, 
a linear combination of the two still satisfies a conserva- 
tion law, namely 
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Similarly, a linear spin-voltage profile is possible when 
there is direct spin relaxation of the conduction electrons, 
but not of the impurity spins. While the spin relaxation 
for the conduction electrons spoils the conservation laws 
for the spin and the spin heat currents, the heat current 
remains conserved. 

In the full problem where direct spin relaxation is 
included for both the impurity spins and the conduc- 
tion electrons, only the charge current is conserved, and 
strictly speaking, thermoelectric generation of a linear 
spin-voltage profile becomes impossible. Nevertheless, it 
is evident from the above discussion that the various cur- 
rents differ in the robustness of their conservation laws, 
and consequently some of the length scales for the expo- 
nential decays can significantly exceed the spin diffusion 
length. As a result, one may still find quasilinear spin- 
voltage profiles in practice even when a linear profile is no 
longer possible. Specifically, we find from a lengthy but 
straight-forward analysis that the three currents J sp i n , 
Jq, and JQ sp in involve the length scales 
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with •S'i = (Sf±Si). Eq. (fT^j) holds for low temperatures, 
ksT <C €p. and weak ferromagnets, /i, <C ey. 
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The length £\ is the shortest of the three length scales 
and corresponds to the usual spin diffusion length. The 
length scale £2 is of the same order as the spin diffu- 
sion length, except when direct spin relaxation is weak 
compared to spin-flip scattering on magnetic impurities. 
In the latter situation, £2 becomes much larger than the 
spin diffusion length and diverges in the absence of di- 
rect spin relaxation. Importantly, £3 is generically much 
larger than the spin diffusion length for weak ferromag- 
nets. According to Eq. (|T2|). the ratio £3/^1 is of or- 
der (ef/e^iT^/Tf) 1 / 2 > 1 when direct spin relaxation 
dominates and (e//e i )(r^/r^ ir ) 1 / 2 > 1 in the opposite 
limit when direct spin relaxation is weak. The length 
scale £3 diverges when no spin relaxation is present for 
electrons or for impurity spins, as well as in the absence 
of inelastic spin-flip scattering. 

Spin voltage. — The largeness of £3 allows for the ther- 
moelectric generation of a quasilinear spin voltage pro- 
file. To illustrate this result and to assess the role of 
the spin thermopowcr S-, we consider a fcrromagnet of 
length L with a constant temperature gradient applied 
along the x-direction, similar to the experimental setup 
in |4J. We then determine the spin- voltage profile from 
Eqs. (fTJ) and (0 together with the boundary conditions 
J S p in O = ±L/2) = 0, T+{x = ±L/2) = ±AT/2, and 
T-(x = ±L/2) = 0. 

Without direct spin relaxation, £2 diverges in addition 
to £3. For low temperatures fc^T -C £f and weak fer- 
romagnets h, <C e/, we then find the spin chemical 
potential profile 



p- = -e (S-T+ + S+T_) 



(13) 



in the bulk of the sample with T+ ~ AT(x/L) and T_ ~ 
T + [(eS + /2~f)e l Tf-2T s l-e l Tf'}/2T s l. Note that both S- 
and S+ contribute to the spin chemical potential profile, 
and it will in general depend on the specifics of the system 
which of the two terms dominates. 

In the general case, with nonzero direct spin relaxation 
for both the conduction electrons and the impurities, we 
can focus on the quasilinear contribution decaying on the 
largest length scale £3 and obtain 
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with T+ ~ ATsinh(a;/£3)/2sinh(L/2£3) and T_ ~ 
T + [(e5 , + /2 7 )e i r+ - 2r s _ f - eiTf'}/2(Tf + Tf T ). Interest- 
ingly, there is no longer a contribution which involves S- . 
The reason for this is that the contribution of S— is sub- 
dominant in the sense that its spatial profile is controlled 
by £2 rather than £3. We therefore find that generically, 
the spin thermopower does not contribute significantly 
to the spin Seebeck effect. 

Conclusions. — We have shown that a (quasi)linear 
spin-voltage profile can be generated thcrmoclcctrically 
in a fcrromagnet with paramagnetic impurities. In the 



absence of direct spin relaxation of conduction electrons 
or impurities, the spin chemical potential acquires a 
strictly linear component due to conservation of a com- 
posite current made up from the heat and spin currents. 
More generally, this conservation law is only approximate 
and a quasilinear behavior of the spin chemical potential 
emerges. Surprisingly, the resulting (quasi)linear spin 
voltage is directly related to the spin thermopower only 
when neither impurities nor conduction electrons are sub- 
ject to direct spin relaxation. While it is not clear to 
which degree our model results apply directly to the re- 
cent experiments, they show that a full understanding of 
the observed quasilinear dependence of the spin chemical 
potential may also involve purely electronic mechanisms. 

The underlying physics of the quasilinear dependence 
of the spin chemical potential is the spin asymmetry 
of inelastic spin-flip scattering caused by partially spin- 
polarizcd magnetic impurities. Both the spin asymmetry 
and the inelasticity should also be a feature of scattering 
of electrons on magnons in itinerant ferromagnets. For 
this reason, our results suggest that it would be interest- 
ing to extend the kinetic theory to include the effects of 
magnons which however, unlike impurities, have a finite 
dispersion and propagate. 

Our analysis also neglects electron-electron interac- 
tions. On the one hand, these will modify the response 
functions <7 CT! and S a , e.g., by intoducing off-diagonal 
spin-drag components. In addition, they provide another 
decay channel for the spin heat current which decays on 
the shortest length scale even in the absence of electron- 
electron scattering. We therefore expect that electron- 
electron interactions modify our results only in quantita- 
tive detail. We hope to return to these issues in a future 
publication. 
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